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ABSTRACT 

We present the line of sight radial velocity dispersion profile of the Sextans dwarf spheroidal 
galaxy (dSph), based on a sample of 88 stars extending to about 1° (~ 1.5kpc). Like the Draco 
and Ursa Minor dSphs, Sextans shows some evidence of a fall-off in the velocity dispersion at 
large projected radii, with significance p - 0.96. Surprisingly, the dispersion at the very centre 
of Sextans is close to zero (with significance p - 0.96). We present evidence which suggests 
that this latter change in the stellar kinematics coincides with changes in the stellar popula- 
tions within the dSph. We discuss possible scenarios which could lead to a kinematically and 
photometrically distinct population at the centre of Sextans. 

Key words: dark matter — galaxies: individual (Sextans dSph) — galaxies: kinematics and 
dynamics — Local Group — stellar dynamics 



1 INTRODUCTION 

It has long been known that the velocity dispersions of Local Group 
dwarf spheroidals (dSphs) imply mass-to-light ratios M/L of up to 
~ 100 Mq/Lq. Until recently, all estimates of dSph M/L have been 
based upon a measurement of the central velocity dispersion, and 
the assumption of mass follows Hght. The advent of datasets of hun- 
dreds of discrete radial velocities of bright stars in the nearby dSphs 
has changed all this. To date, the line-of-sight velocity dispersion 
profiles in the Fornax, Draco and Ur sa Minor dSphs have been 
mapped out right to the op tical edge jMatedll99i iKlevna et all 
l2002tlWilkinson et alllOM) . 

In this paper we turn our attention to the Sextans dSph and 
present the first results from our new radial velocity survey. The 
new data show an almost flat velocity dispersion to large radii, with 
some evidence of a fall- off in the outer parts si milar to that seen in 
Draco and Ursa Minor (" Wilkinson etd]l2004 . Surprisingly, Sex- 
tans appears to be kinematically cold at the very cent re. This is the 
secon d localised cold population observed in a dSph iKlevna et alj 
l2002h . and it coincides with observable changes in the composition 
of the stellar population. We suggest that the nucleus of the Sextans 
dSph contains a dynamically separate component, possibly the re- 
mains of a star cluster dragged to the centre by dynamical friction. 



2 OBSERVATIONS 
2.1 Velocities 

We observed Sextans with the multifibre instrument AF2/WYFF0S 
on the William Herschel Telescope on La Palma on 11-16 March 
2004. Approximately two of these six nights were lost to weather 



or technical problems. We reduced the spectra using a modified 
version of the IRAF WYFRED task, and we obtained the final ve- 
locities using the FXCOR cross-correlation function (CCF) veloc- 
ity package. From the point of view of the data reduction. Sex- 
tans presents particularly thorny difficulties. Due to the bulk line of 
sight velocity of Sextans, all three Ca triplet lines absorption lines 
used to measure the velocity are embedded in bright sky emission 
lines. Hence it is of paramount importance to subtract the sky spec- 
trum accurately, and to distinguish between those correlation peaks 
caused by stellar absorption lines and those arising from a random 
or systematic sky residual. 

Each WYFFOS pointing consisted of approximately 75 fibres 
allocated to objects, and an average of 27 sky fibres interleaved 
among the object fibres. After rebinning each spectrum to a com- 
mon linear dispersion and flattening with a tungsten lamp spectrum, 
we constructed a sky spectrum for each object fibre by median- 
combining the five sky spectra closest to the object spectrum on 
the CCD, thereby minimising the effects of dispersion-direction 
PSF variations. Next, we subtracted the sky spectra from the ob- 
ject spectra by modelling each spectrum as the sum of a third-order 
polynomial plus a multiple of the sky, and minimising the inte- 
grated absolute deviation of the model from the actual spectrum, 
excluding those regions known to contain stellar absorption lines. 
Effectively, this procedure subtracts out the correct amount of sky 
by minimising the sky line residuals. As the final step before com- 
puting velocities, we shifted each spectrum to a common heliocen- 
tric velocity frame. 

For our previous work on Dra co and UMi ( Klevna et all200l|. 
l2002ll2003HWilkinson et alJl2004) . the sky lines were safely sep- 
arated from the stellar lines, and it was possible to combine all of 
the good spectra for an object into a single spectrum, producing an 
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unambiguous CCF against a synthetic template. With Sextans, it 
was difficult to distinguish good spectra from bad, so we elected to 
obtain as many individual measurements as possible by measuring 
a velocity with each of the three Ca triplet lines in each one hour 
exposure, using FXCOR in an automated mode, and allowing it to 
select the highest CCF peak between -200 and +400 kms"'. For 
each velocity, FXCOR also giv es an uncertainty oc ( 1 where 
R is the Tonry-Davis R value ( Tonrv & Davis 1979). 

Next, we constructed the following sigma-clipping procedure: 
for a set of measured velocities and their nominal uncertainties cr,, 
we compute the median of the set, discard the farthest outlier (in 
terms of cr,), and iterate until all of the remaining points are less 
than 2.5 cr, from the median. Then we computed the velocity of 
each object in each of lines 1, 2, and 3 by applying the sigma clip- 
ping procedure to all measurements in the line in question. We ob- 
served that the velocities in line 2 and 3 were consistent, but gen- 
erally disagreed with line 1 - this is not surprising, as line 1 is the 
weakest line, and is atop the strongest sky line. Thus we elected 
to compute our final best velocity for each object by applying the 
sigma clipping procedure to the combined set of velocities from 
lines 2 and 3. 

To accept a velocity for our subsequent analysis, we required 
that there exist at least four velocities, representing at least 2/3 of 
the inputs, which survive the sigma clipping procedure. In all, 88 
Sextans member stars from the combined 2 + 3 lines within 30 
km s"' of the dSph's median velocity met these criteria. 

Our sigma clipping and subsequent error estimates depended 
on the validity of the i?-value velocity uncertainties produced by 
FXCOR. Each final velocity, as a Gaussian-weighted average of in- 
dividual velocities, had an associated^-. The mean;^'^ per degree of 
freedom for all member and non-member stars with good velocities 
was 0.93, suggesting that the FXCOR velocities were indeed valid. 
Repeating the clipping procedure after rescaling the FXCOR errors 
by 0.75 and 1.33 produced dix^ per degree of freedom that deviated 
from the expected value of 1 in the expected directions. The small- 
est Pix^) observed for a Sextans member was 0.004, which would 
be expected to occur by chance in a sample of this size a third of 
the time. 

Despite our care to sift out spurious velocities from our data, 
the sky line contamination remains a concern. Fig.Qshows the re- 
sult of reducing our sky fibres as though they were object fibres. 
Because lines 1 and 2 of the Ca triplet are directly under two strong 
sky lines, the residue of the sky tends to produce spurious veloc- 
ities at 233 and 223 kms"', respectively. The peaks actually be- 
comes stronger if we limit ourselves to larger Tonry-Davis R val- 
ues, probably because we are picking out spectra with systematic 
sky subtraction problems. Only line 3, embedded in a dense cluster 
of lower amplitude sky lines, shows an unbiased distribution of spu- 
rious velocities. Because of this potential source of contamination, 
we will take care to verify any observed kinematical qualities of 
Sextans using line 3 alone, as well as with our best 2 + 3 velocities. 

Fig. El shows our individual velocities as a function of pro- 
jected radial distance. The data using only line 1 are obviously of 
much poorer quality than for lines 2 and 3; in particular, the median 
velocity for line 1 is 232.5 kms" , whereas the medians for lines 
2 and 3 are 226.6 and 225.3 km s"', respectively. Immediately ap- 
parent to the naked eye is the apparent tendency of the outermost 
stars of the top 2 + 3 panel to cluster very tightly around Sextans' 
mean velocity; this kinematic coldness at large R was also fou nd to 
be present in UMi and possibly Draco jWilkinson et all2004h . Un- 
fortunately, the data for lines 2 and 3 alone are not of sufficiently 
good quality to support this conclusion on their own. 
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Figure 1. Distribution of pseudo-velocities of sky fibres when processed as 
though they were data fibres, for each of the three lines of the Ca triplet. 
Lines 1 and 2 produce a spurious peak near the mean velocity of Sextans, 
resulting from random and possibly systematic sky residuals. Line 3, on the 
other hand, produces a relatively unbiased distribution of spurious veloci- 
ties. 



2.2 Comparison with previous work 

In the past, repeated velocity measurements of individual Sextans 
stars have been inconsistent from observer to observer, probably 
because of the severe sky contamination. Suntzeff et al. (1993), 
(henceforth S93), measured 43 velocities in Sextans using the Ar- 
gus multifibre spectrograph on the CTIO 4m; five of these were in 
the previous set of Da Costa et al. (1991). On the basis of the dis- 
agreement of the two data sets, S93 argue that Da Costa et al. (1991) 
underestimated their errors. Similarly, Hargreaves et al. (1994), 
hence H94, observed 26 Sextans spectra using the single-slit ISIS 
instrument on the WHT; 1 1 of these were in the sample of Suntzeff 
et al. (1993). H94 find that the two data sets disagree, and conclude 
that S93 understated their velocity uncertainties. 

Continuing the Sextans tradition of finding fault with our pre- 
decessors, we find that our data agree neither with S93 nor with 
H94, with whom we share 16 and 7 stars, respectively. For the pur- 
pose of comparing our velocities to previous results, we consider 
only the 38 S93 velocities without a "very poor correlation", and 
the 21 H94 velocities that have a Tonry-Davis R above their cutoff 
of 7.5. We find a discrepancy, expressed in^^ per degree of free- 
dom, ofx^/Niof = 2.44 with S93, comparable to the;(^V-^Vdof - 2.82 
observed between S93 and H94. With H94, we find a disagree- 
ment of ;t'^/A'dof = 7.46; this value is larger than the disagreement 
with S93 probably because the errors of H94 are much smaller 
(2 km s"') than those of S93 (median 5.5 kms"'). 

Figure|3|shows all stars that overlap in at least two of the data 
sets. None of the three data sets is in obvious systematic disagree- 
ment with the other two. Compared with S93, we have the advan- 
tage of a more stable instrument and a more robust wavelength so- 
lution - S93 needed to adjust their final velocities by a large system- 
atic factor dependent on the fibre number, whereas our dispersion 
solution was stable to ~ 0.8 kms"', verified by cross-correlating 
each non sky-subtracted spectrum against a combined sky template. 

Compared with H94, we have the benefit of many more mea- 
surements per star, which we increased by a further factor of three 
by using each line separately; this permits testing each star's in- 
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Figure 2. Individual velocities as a function of projected radius R for the 
combined data (using lines 2 and 3 of the Ca triplet), and for lines 1, 2, 
and 3 individually. Outer dotted lines show our membership limits, and the 
inner dotted line is the median of the entire data set. 



ternal consistency. H94, on the other hand, occasionally use only 
one measurement per star, relying on a global calibration of error 
as a function of Tonry-Davis R. On the other hand, we accept low 
values of Tonry-Davis R, opting to reject outliers with the proce- 
dure described in ^2.11 As described above, this procedure gives a 
final Pix^) very consistent with correct normally distributed Tonry- 
Davis R errors. In further Monte-Carlo tests of Sextans-like data 
50% contaminated with a flat distribution of spurious velocities, it 
has proven to be very robust at removing the outliers and produc- 
ing data sets with the correct mean and dispersion. Nevertheless, 
concern remains that both we and previous observers might be per- 
mitting spurious or sky-shifted correlation peaks into the final data 
set. 



2.3 Velocity dispersion 

For their entire samples, S93 and H94 reported velocity disper- 
sions of 6.2;^q ' kms" 



and 7.0^[ q kms"', respectively. However, 
using their more appropriate maximum likelihood formulation for 
the dispersion rather than the method of Armandroff & Da Costa 
(1986), H94 recompute the dispersion of the S93 data to be between 
6 and 8 kms"', depending on the cut chosen for the Tonry-Davis 
R value. 

H94 express the probability of the line-of-sight velocity dis- 
persion (v^) and the mean velocity v for a sample of A' velocities v, 
with errors ct; as 



N exp 



1 (Vj-v)- 



(1) 



They then find the optimal value of (v") and v iteratively, and adjust 
(v^) by a factor of A'(A' - 1)"' because the number of degrees of 
freedom is reduced by one through the simultaneous computation 
of V. Finally, they expand P((v->''^, v) around the optimal point in 
(i'^>, V to obtain the errors. 

We note, however, that Eq.Qis a Gaussian in v, so it can be in- 
tegrated analytically over v to obtain a one dimensional probability 
P«v2)'/2): 
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Figure 3. Subset of Sextans velocities (with l-cr en'or bars) common to at 
least two of the following: this work, Suntzeff et al. (1993), or Hargreaves 
et al. (1994). Each star occupies a unique position along the horizontal axis, 
with the two or three different measurements along the vertical axis. 



2;r exp[-|(c-fc2/a)] 

^ r 

« nr27r(ff-2 + <v2))"' 



(2) 



where a = 2;«v^) + crf)-\ b = 2, v,«v2) + crj)-\ and c = 
2, v^((v^) + cr^)"'. P({v^)) can then be integrated numerically to 
obtain precise error bounds. This is the approach taken in this pa- 
per: when collectively computing the dispersion of an entire data 
set, we use the P({v^)) obtained from Eq.|2| effectively integrating 
over all values of v. However, when binning the data radially and 
comparing dispersions among bins, we wish to use a single mean 
velocity across all bins, so we instead compute P({v^)\v) by fixing 
V in Eq. Q at the median velocity of our entire data set. In these 
cases, we may introduce the further modification that the Gaussian 
in Eq.Qis replaced by a Ga ussian convolved wi th a binary velocity 
distribution, as described in 'Klevna et aljj2002h . To ensure that our 
implementations of Eq.Qand Eq.|2|are correct, we have conducted 
extensive Monte-Carlo tests using artificial data. 

We note that the data of S93 are contained entirely within the 
central 16.35', whereas the velocities of H94 extend outside this 
radius. For consistency, we limit our comparisons among different 
observers' dispersions to this common central region. Using our 
global, purely Gaussian P{{v^)), we compute the following most 
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Figure 4. Velocity dispersions inside bins with boundaries at projected ra- 
dius R = 5', 10', 15', 20', 25', 30', and 40'. The top panel has velocities 
computed using both Ca triplet lines 2 and 3, and the other panels are for 
lines 1 , 2, and 3 individually. Velocities using line 1 are of poor quality 
and give much larger dispersions than vel ocities obtained wit h lines 2 and 
3. A binary fraction of 40% is assumed jKlevna et alJl20o3) . Each point 
represents the most likely dispersion, and the error bar is the central 68% 
confidence interval. The R value of each point is the mean projected radius 
of the stars in the corresponding bin. The three good data sets are suggestive 
of a dispersion that is zero in the centre, rises to about 8 kms"', and then 
falls off at large radii. 



likely dispersions inside 16.35', with central 68% confidence in- 



tervals - this paper: (v ) 
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7.80+ 



km s 



from 38 stars; H94: 8.62+f^3 kms" 



' from 30 stars; S93: 
from 14 stars. 

These dispersions are somewhat larger than the 6 to 7 kms"' 
dispersions given by S93 and H94. Nevertheless, they are statis- 
tically consistent, and the disparities probably arise from the differ- 
ent subsamples used, and the slightly different (and perhaps more 
rigorous) method we use to compute the dispersion. We also ex- 
clude the low-dispersion stars of H94 outside 16.35', increasing the 
dispersion we measure in the H94 data. Given the disagreements 
we noted among the three data sets, it is reassuring that they give 
such similar dispersions. This agreement may arise from the fact 
that the dispersion is sensitive to excess velocity errors only when 
these errors become large compared to the dSph's overall disper- 
sion, whereas the disagreement among data sets becomes large 
when the excess velocity errors are large compared to the nominal 
velocity uncertainties. 

Fig-Elshows the final velocity dispersions as a function of pro- 
jected radius for our data set. From Fig.|4| it is difficult to evaluate 
whether the statistical significance of the apparent fall in the dis- 
persion at small and large R is an artefact of the adopted binning. 
Accordingly, in Fig. |5|we plot the velocity dispersion advancing 
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Figure 5. Velocity dispersion computed using the N = 1 ... 40 innermost or 
outermost stars, for our combined 2 + 3 dataset. The point at N represents 
the most likely dispersion for the sample of stars from I to N using the 
purely Gaussian formulation of Eq.Qand setting v equal to the median of 
our entire sample. The bars are the central 68% confidence intervals at each 
point. The bars need not overlap the most likely point if the probability 
profile is asymmetric. 
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Figure 6. Probability that the dispersion in the N innermost or outermost 
stars is larger than the dispersion of the entire 88 star 2+3 sample, computed 
according to Eq.|3 There is > 2o" evidence of a dispersion drop in both the 
centre and outer limits of Sextans. 



from the outermost or innermost star, each time adding one star 
to the sample. We find that the most likely dispersion for the four 
outermost and seven innermost stars is zero. As we move inward 
or outward adding stars to the sample, the upper bound on the dis- 
persion falls until we reach a region of higher dispersion. Next, we 
compute the probability that the inner and outer samples have a 
dispersion smaller than the entire ensemble of Sextans stars as fol- 
lows: if we have two samples of data a and b with a common mean 
V, the joint probability distribution for their dispersions is 

P{{^HJ\ {V')T) = J P«y-)a\ V) X Pav')l'\ V) d V (3) 

where P{{v^)'J^ ,v) is given in Eq.Q The probability that {v^)'J^ > 

9 1/2 

(v )j is then 
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Figure 7. V,V ~ I colour-magnitude diagram of stars within 20 arcmin of 
the centre of Sextans (based on archival INT images). Over-plotted are the 
5 bins used to define the stellar sub-populations. See text for a detailed 
discussion. 



This is accomplished identically to the step that converts Eq.Qto 
Eq. 121 Fig. |6| shows the probability, using Eq.|4| that the A' inner- 
most or outermost stars have a velocity dispersion smaller than that 
of the rest of the sample. We find that the 17 outermost stars of 
Sextans are kinematically colder than the rest of the dSph at up 
to the p = 0.99 confidence level. We can draw the cutoff almost 
anywhere between 4 and 17 and obtain the same result, so this con- 
clusion is not finely tuned to the value of our one free parameter, 
the cutoff radius. Similarly, the innermost 7 stars are also cold at 
the p = 0.99 confidence level. In this case, the choice of cutoffs 
is narrower, with only 4 to 7 giving this level of confidence. We 
also performed a Monte-Carlo test using simulated data with a flat 
8 km s"' dispersion and our nominal observational errors to deter- 
mine how frequently we would expect to observe dips of the mag- 
nitude seen in Fig. |6| anywhere in the first or last N = 24 stars - 
we find that the lowest outer and inner dispersion dips are signifi- 
cant at the p = 0.96 level. Recreating Fig.|6|for Ca triplet lines 2 
and 3 taken alone gives similar dips, but at lower significance levels 
p = 0.8 to p = 0.95. 



3 MODELLING 

One possible explanation for the observed data would be the pres- 
ence of a kinematically distinct stellar component at the centre of 
Sextans. To test this idea further, we use the procedure described 
m lllevna et al. 1 2003) to examine the position- velocity data set for 
signs of substructure. We scan the face of the dSph with an aperture 
of radius 5'. At each aperture location, we calculate the likelihood 
that these velocities are drawn from (i) a single Gaussian with dis- 
persion (v^) = 7.97 kms"' (ii) a velocity distribution in which a 
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Figure 8. Top: Suii'ace brightness profile for Sextans together with the best- 
fitting single-component Plummer law (lower curve) and a Plummer law 
with a central power-law cusp of index 4.4 (upper curve). The two curves 
deviate only inside about 1.6 arcmin where the cusped model rises to fit 
the innermost data point. Observed data taken from llrwin & HatzidimitrioiTI 
O995). Bottom: Velocity dispersion profile with best-fitting model (cenfi'al 
data point omitted from the fit). 
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Figure 9. Background-subtracted projected density profiles of the red (solid 
curve) and blue (broken curve) horizontal branch stars (top) and of the red 
(solid curve) and blue (broken curve) RGB stars (middle) in Sextans based 
on INT photometric data. The bottom panel shows the radial variation of 
the HB morphology index. All three panels show that the composition of 
the stellar population at the centre of Sextans varies with radius on a scale 
of about 10 arcmin. See text for a detailed discussion. 



fraction / of the stars are drawn from the main Gaussian and the 
remainder belong to a kinematically distinct population with differ- 
ent mean velocity and dispersion. Only in an aperture centred close 
to the origin does any model other than the main Gaussian represent 
a significantly better fit to the velocity data. A model in which all 
the stars in this aperture belong to a population with velocity dis- 
persion 0.5 km s" (and a mean velocity within 2 km s" of the bulk 
motion of the dSph) is about 195 times more likely than a model 
in which they belong to the main population. Monte Carlo realisa- 
tions of our data set demonstrate that the likelihood of observing 
a spurious cold clump with the same probability ratio at any loca- 
tion in the cluster is relatively high. However, only 4.6 percent of 
realisations display a clump with an offset in mean velocity of less 
than 5 km s"' . We therefore conclude that the central region shows 
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evidence for the presence of a kinematically distinct population at 
about the 2cr confidence level. 

The light distribution of Sextans displays a sharp central 
rise (see| lrwin & Hatzidi mitrio u 1995, and the top panel of Fig.|8}. 
It is tempting to imagine an association between this feature and the 
cold population. However, the difference in the length scales of the 
two features means that a direct association is very unlikely. At a 
radius of 2 arcmin (the radius of our innermost velocity data point) 
the ratio of the central component to the main component is already 
less than unity, for any plausible power-law fit to the centre of the 
light distribution. Thus, if the central density profile were represen- 
tative of the spatial distribution of the cold population, we would 
not expect that six of the innermost eight stars in our velocity sam- 
ple would be cold. 

A number of authors have presented evidenc e for the pres- 
ence of radial stellar population gradients in Sextans (Harbeck et al. 
[2001 ; Lee et al. 2003). In particular, the red horizontal branch 
(RHB) stars appear to be more centrally concentrated than the blue 
horizontal branch (BHB) stars. Following these authors, we define 
five bins in the (V, V-I) colour magnitude diagram (CMD) obtained 
from deep {V < 22.5) INT imaging of Sextans (see Fig.0. There 
are three bins on the horizontal branch (RHB, BHB, RRLyrae) 
and two bins on the red giant branch (RGB) between V=18.8 and 
V=19.8. The red RGB bin (RRGB) traces the mean colour of the 
RGB while the blue bin (BRGB) covers the blue side of the RGB. 
The number of stars in each bin is corrected for foreground con- 
tamination by counting the number of stars in a bin of equal area 
displaced away from the CMD features associated with Sextans. 

In the top panels of Figure |9| we compare the background- 
subtracted surface density profiles of the red and blue compo - 
nen ts of th e HB and RGB. In agreement with Harbeck et al. ( 2001) 
and lLee et aL 1.2003.) , we find that the RHB is significantly more 
centrally concentrated than the BHB. There is also weak evidence 
that the RRGB is more concentrated than the BRGB. The bottom 
panel of the Figure shows the radial variation of the HB morphol- 
ogy index (;ibhb - "rhb)/(«bhb + "rhb + "rrl), where ;ibhb, «rhb 
and wr rl are the numbers of BHB, RH B and RRLyrae stars, respec- 
tively iLee. Demaraue. & Zinnll990ll . Within 10 arcmin, the HB is 
dominated by its red component, while at larger radii, the numbers 
of red and blue HB stars are comparable. Given the uncertainties in 
the definitions of the CMD bins in our analysis and the unknown 
level of contamination of the RHB bins by RRLyraes, the details 
of these radial variations are not well determined. However, it is 
clear that the composition of the stellar population at the centre of 
Sextans varies on a scale of about 5-10 arcmin. 

If both the stellar components of Sextans are tracer popula- 
tions in a dark matter dominated potential well, then the more con- 
centrated population should have a smaller velocity dispersion. The 
RGBs of the two populations may be coincident in the CMD and 
so we would not necessarily expect to see correlations in our ve- 
locity data with position on the RGB. However, we note that in the 
innermost radial bin, the surface density of RHB stars is a factor 
of 9 higher than that of the BHB stars. It is therefore plausible that 
within 5 arcmin of the centre we would observe mostly stars asso- 
ciated with the cold component. 

One possible origin for the cold, inner population of Sextans 
is the preferential retention of star-forming gas at the centre of the 
dSph leading to a yo unger and/o r more metal-rich, centrally con- 
centrated population faarbeck e t al. 2001). An alternative explana- 
tion is that the cold central population originated in a star cluster 
which spiralled to the centre of Sextans under the influence of dy- 
namical friction. During the final stages of this process, the clus- 



ter deposited stars in the inner regions of the dSph. It is possible 
that the central cusp in the light distribution is the final, probably 
unbound, remnant of this cluster. A similar scenario has been pro- 
posed by Oh & Lin 1 2000) for the origin of nucleated dwarf galax- 
ies. 

In addition to providing an explanation for the photometric 
and kinematic features discussed already, the cluster scenario nat- 
urally accounts for the distribution of blue stragglers in Sextans. 
Lee et al. (2003) noted that the brighter (i.e. more massive) blue 
stragglers are more centrally concentrated than the fainter ones. 
As those authors comment, mass-dependent spatial distributions of 
stars are a generic outcome of mass segregation in stellar systems. 
However, the time scale for such segregation to occur in Sextans 
is very long due to the low surface brightness of the stellar popu- 
lation. If a significant fraction of the blue stragglers were formed 
in a star cluster which subsequently disrupted near the centre of 
Sextans, mass segregation within the cluster would ensure that the 
most massive blue stragglers would be the last to be tidally removed 
from the cluster and hence would have a more concentrated spatial 
distribution. 

If we assume that all the light in the central 10 arcmin is due 
to cluster debris (as suggested by the dominance of the cold popu- 
lation in the inner regions and the large ratio of RHB to BHB stars 
inside 10 arcmin) and that the (con stant) luminosity density within 
10 arcmin is 0.002Lopc-3 lMatMl ll998h . we obtain a luminosity 
of 1.3 X 10^ Lq for the putative cluster. If the mass to light ratio is 
about 2 solar units (typical of Galactic globular clusters), then the 
mass of the cluster is about 3 x IO^Lq. This is comparable to the 
masses of Galactic globular clusters. 

We determine the mass profile of Sextans in the context of 
the cluster model by assuming a Plummer model for the light dis- 
tribution, excluding the central point (see Figure|HJ. Similarly, we 
remove the central point from the velocity dispersion profile and fit 
a profile of the form (see Figure|H) 



<v'>o 



(1 -I- (R/R^yyi'' 



(5) 



If we further assume isotropy of the velocity dispersion tensor, 
then the mass within 1 kpc (using the Jeans equation) lies in the 
range 3 x 10^ - 1.5 x IO^'Mq ~ the exact value is sensitive to the 
form assumed for the outer fall-off in the dispersion. Assuming 
that the dark matter has an isothermal distribution out to 1 kpc, 
then the time taken for a cluster of mass 3 x IO^'Mq to spiral to 
the centre due to dynamical friction would be in the range 0.7-1.5 
Gyr I Binnev & Tremaine 1987). Thus, we would expect to find any 
such cluster at the centre of the Sextans dSph. 

We also note that the mass within 10 arcmin (= 250 pc) 
in this model is approximately 2.5 x IO^Mq. Hence, the im- 
plied mass to light ratio of the inner regions is roughly 19 so- 
lar units. This is much higher than the M/L expected for an 
old stell ar population based on observations of globular clus- 
ters (e.g. lParmentier & Gilmorj|200lLlFeltzing. Gilmore. & Wvsd 
1999}). This suggests that the inner regions of Sextans are dark mat- 
ter dominated, which is consistent with our earlier assumption that 
the stellar populations in the central regions are not currently self- 
gravitating. 



4 CONCLUSIONS 

We have presented the radial variation of the velocity dispersion 
for the Sextans dSph. The behaviour of the velocity dispersion pro- 
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file at large radii is qualitative ly similar to that seen in the Draco 
and Ursa Minor dSphs I Wilki nson et alj Eo04y In those galaxies, 
the fall-off in the velocity dispersion at large radii coincided with a 
break in the light distribution and was interpreted as indicative of 
mild tidal perturbation of the stars. In Sex tans, the photometry pre- 
sented bv llrwin & Hatzidimitriou | il995h is of insufficient quality 
to identify features at 40 - 50 arcmin. The interpretation of the fall- 
off in the outer velocity dispersion of Sextans is unclear, but may 
be due to the tidal field of the Milky Way. 

Intriguingly, the dispersion at the very centre of Sextans is 
close to zero and this low value is coincident with significant radial 
gradients in the stellar populations. We suggest that this is caused 
by the sinking and gradual dissolution of a star cluster at the cen- 
tre of Sextans. This is consistent with both our inferred mass for 
Sextans and the estimated timescale for dynamical friction to bring 
any star cluster to the centre of the Sextans. The model also pro- 
vides an explanation for the origin of the central cusp in the light 
distribution and the different spatial distributions of the bright and 
faint blue stragglers in Sextans. Thus, our hypothesis provides a 
reasonable explanation of all the observed data, both kinematic and 
photometric, on the Sextans dSph. 



Tonry J., Davis M. 1979, AJ, 84, 151 1 

Wilkinson M.I., Kleyna, J.T., Evans, N.W., Gilmore G. 2001, MN- 
RAS, 330, 778 

Wilkinson M.I.,Kleyna J.T, Evans N.W., Gilmore G., 2004, ApJ, 
999, L999 



ACKNOWLEDGEMENTS 

JTK gratefully acknowledges the support of the Beatrice Watson 
Parrent fellowship. MIW acknowledges the support of PPARC. We 
are very grateful to Mike Irwin for providing the astrometric cat- 
alogue used to produce our target list. We thank George Seabroke 
and the staff" at the Isaac Newton Group at La Palma for their assis- 
tance in obtaining the data for this paper and Dougal Mackey for 
useful discussions. We are grateful to the referee Carlton Pryor for 
useful suggestions. 



REFERENCES 

Armandroff T.E., Da Costa G.S. 1986, AJ, 92, 777 

Binney J., Tremaine S. 1987, Galactic Dynamics, Princeton Uni- 
versity Press, Princeton 

Da Costa G. S., Hatzidimitriou D., Irwin M. J., McMahon R. G. 
1991, MNRAS, 249, 473 

Feltzing, S., Gilmore, G., & Wyse, R. R G. 1999, ApJ, 516, L17 

Harbeck, D., et al. 2001, AJ, 122, 3092 

Hargreaves J. C, Gilmore G., Irwin M. J., Carter D. 1994, MN- 
RAS, 269, 957 [H94] 
Irwin M., Hatzidmitriou D., 1995, MNRAS, 277, 1354 
Kleyna J.T, Wilkinson M.I., Evans N.W., Gilmore G. 2001, ApJ, 
563, LI 15 

Kleyna J.T, Wilkinson M.I., Evans, N.W., Gilmore G. 2002, MN- 
RAS, 330, 792 

Kleyna J.T, Wilkinson M.I., Gilmore G., Evans N.W. 2003, ApJ, 
588, L21 

Lee, M. G., et al. 2003, AJ, 126, 2840 

Lee, Y., Demarque, R, & Zinn, R. 1990, ApJ, 350, 155 

Mateo, M., 1997, in The Nature of Elliptical Galaxies; 2nd 

Stromlo Symposium, ASP Conf. Ser. 116, 259 
Mateo, M. L. 1998, ARA&A, 36, 435 
Oh K.S., Lin D.N.C 2000, ApJ, 543, 620 
Parmentier, G. & Gilmore, G. 2001, A&A, 378, 97 
Suntzefi" N.B., Mateo M., Temdrup D.M., Olszewski E.W., 

GeislerD., Weller W. 1993, ApJ, 418, 208 [S93] 



